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Abstract 


A  new  turbulence  model  is  used  to  describe  the  acoustical  scattering  from 
atmospheric  turbulence.  A  complete  set  of  fluid  equations,  including  the  heat 
flow  equation  with  zero  conductivity,  is  presented  for  an  ideal  gas  atmosphere. 
From  this  set,  a  complete  set  of  coupled  linear  diff'erential  equations  is  de¬ 
rived  for  the  acoustic  pressure,  temperature,  mass  density,  and  velocity  in 
the  presence  of  stationary  turbulence.  From  these  acoustic  wave  equations, 
expressions  for  acoustic  scattering  cross  sections  are  derived  for  individual 
localized  stationary  scalable  turbules  of  arbitrary  morphology  and  orienta¬ 
tion.  Averages  over  random  turbule  orientations  are  also  derived.  Criteria  for 
comparability  of  orientationally  averaged  turbules  with  different  envelope 
functions  are  presented  and  applied,  and  cross  sections  for  Gaussian  and  ex¬ 
ponential  envelopes  are  compared.  The  azimuthal  dependence  of  the  velocity 
scattering  cross  section  for  a  spherically  symmetric  nonuniformly  rotating 
turbule  is  illustrated.  It  is  shown  that,  for  incoherent  scattering,  a  collection 
of  randomly  oriented  turbules  of  arbitrary  morphology  may  be  replaced  by 
an  “equivalent”  collection  of  spherically  symmetric,  nonuniformly  rotating 
turbules  with  randomly  directed  rotation  axes. 


Contents 


1  Introduction  1 

2  Basic  Theory  4 

2.1  Fluid  Equations .  4 

2.2  Acoustic  Equations .  5 

2.2.1  Linear  Acoustics  .  5 

2.2.2  First-Order  Acoustic  Equations .  6 

2.3  Scattering .  ^ 

2.3.1  General  Formulation .  7 

2.3.2  Born  Approximation .  9 

3  Scattering  from  Model  Turbules  15 

3.1  General  Expressions  .  16 

3.2  Isotropic  Ensembles .  18 

3.2.1  Orientation  Averages .  18 

3.2.2  Comparable  Turbules .  19 

3.2.3  Examples . 20 

3.3  Individual  Turbule  Model  .  26 

4  Summary  and  Discussion  29 

References  81 

Distribution  83 

37 

iii 


Report  Documentation  page 


Figures 


1 


2 


3 


4 


5 


6 


Temperature  scattering  efficiency  versus  scattering  an¬ 

gle  9,  for  6T/T00  =  1-0,  for  a  Gaussian  isotropic  ensemble  of 
turbules  of  scale  length  a,  for  several  size  parameters  ka  = 
27ra/A,  A  =  acoustic  wavelength,  in  Born  approximation  .  .  . 

Velocity  scattering  efficiency  at,(^)/7ra^  versus  scattering  angle 
6,  for  V/Coo  =  1-0,  for  a  Gaussian  isotropic  ensemble  of  tur¬ 
bules  of  scale  length  a,  for  several  size  parameters  ka  =  27ra/ A, 
A  =  acoustic  wavelength,  in  Born  approximation . 

Total  Born  approximation  (temperature,  velocity)  scattering 
efficiency  (ar/7ra^,  a^/rra^)  versus  size  parameter  ka  =  2'KalX, 
A  =  acoustic  wavelength,  a  =  turbule  scale  length,  for  a  Gaus¬ 
sian  isotropic  ensemble,  for  5T /Too  =  u/coo  =  1-0  . 

Temperature  scattering  efficiency  '0x9 /‘no?-  versus  scattering 
angle  9,  for  8T/Toc  =  1-0,  for  an  exponential  isotropic  ensem¬ 
ble  of  turbules  of  scale  length  a.  for  several  size  parameters 
ka  =  27ra/A,  A  =  acoustic  wavelength,  in  Born  approximation 

Velocity  scattering  efficiency  0^9/110?  versus  scattering  angle 
9,  for  v/Coo  =  1-0,  for  an  exponential  isotropic  ensemble  of 
turbules  of  scale  length  a,  for  several  size  parameters  ka  = 
27ra/A,  A  =  acoustic  wavelength,  in  Born  approximation  .  .  . 

Total  Born  approximation  (temperature,  velocity)  scattering 
efficiency  (^r/Tra^,  a^/Tra^)  versus  size  parameter  ka  =  27ra/A, 
A  =  acoustic  wavelength,  a  =  turbule  scale  length,  for  an  ex¬ 
ponential  isotropic  ensemble,  for  6T/T00  =  n/coo  =  1.0  .  .  .  . 


21 


21 


22 


24 


25 


25 


IV 


1.  Introduction 


The  motivation  for  the  work  described  in  this  report  was  to  investigate  the 
utility  of  a  new  model  of  atmospheric  turbulence  when  used  to  predict  scat¬ 
tering  of  acoustical  signals  by  the  turbulent  atmosphere.  Much  of  the  in¬ 
formation  in  this  report  we  provide  in  a  journal  article  [1],  in  which  we 
emphasize  the  aspects  of  the  work  strictly  related  to  acoustical  scattering. 
In  this  report,  however,  we  also  include  aspects  of  the  work  not  emphasized 
in  the  article.  For  example,  as  part  of  the  work,  we  examined  in  detail  the 
fluid  flow  equations  to  establish  the  relative  order  of  magnitude  of  the  vari¬ 
ous  terms  retained  when  the  acoustic  wave  equation  is  derived;  for  the  most 
part,  we  did  not  include  this  analysis  in  the  journal  article.  We  also  derived 
the  properties  of  structure  functions  predicted  by  the  new  turbulence  model; 
we  included  in  the  article  only  a  minor  portion  of  these  findings.  Instead,  the 
structure  function  findings  are  recounted  in  a  companion  technical  report  [2]. 
In  short,  all  our  findings  are  included  in  this  report  and  its  companion  [2], 
but  only  those  parts  closely  connected  to  acoustical  scattering  are  included 
in  the  journal  article  [1]. 

Several  different  approximate  wave  equations  for  acoustic  pressure  or  acoustic 
potential  have  been  derived  to  describe  acoustics  in  a  turbulent  atmosphere 
[3-11].  Several  of  these  equations  are  valid  only  to  first  order  in  station¬ 
ary  turbulent  temperature  variation  and/or  velocity.  These  include  Monin's 
equation  [3]  and  its  short- wavelength  limit  [4],  and  Abdullayev's  and  Osta- 
shev’s  [8]  and  Ostashev's  [10]  generalizations  of  Monin’s  equation  to  include 
water  vapor  and  other  additives  and  nonstationarity  of  turbulence.  These 
first-order  wave  equations  usually  provide  an  adequate  basis  for  describing 
linear  atmospheric  acoustics  because  of  two  general  features  of  turbulent 
atmosphere:  first,  the  ratio  of  turbulent  velocity  to  acoustic  wavespeed  and 
that  of  turbulent  temperature  variation  to  ambient  temperature  are  generally 
very  small  for  the  atmosphere,  and  second,  turbulent  time  scales  are  generally 
much  longer  than  acoustic  time  scales.  Section  2  presents  the  basic  theory 
of  acoustic  wave  equations  and  the  scattering  properties  of  turbulence.  Sec¬ 
tion  3  develops  the  scattering  properties  of  model  velocity  and  temperature 
distributions. 
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Section  2  of  this  report  has  four  purposes.  First,  we  provide  a  complete  set  of 
fluid  equations  in  terms  of  pressure,  density,  temperature,  and  fluid  velocity, 
and  a  complete  set  of  coupled  linear  differential  equations  that  describe  linear 
acoustics  in  the  presence  of  turbulence  in  an  ideal  gas  atmosphere.  Second, 
we  show  that  Monin’s  first-order  acoustic  wave  equation  results  from  this 
set.  Third,  we  provide  a  coherent  development  of  general  properties  of  the 
acoustic  scattering  predicted  by  these  wave  equations,  properties  that  are 
independent  of  the  detailed  morphology  of  the  turbulent  velocity,  density, 
pressure,  and  temperature  fields,  but  dependent  on  the  wavelength  and  the 
length  scales  of  the  turbulence  and  the  scattering  volume.  Fourth,  we  estab¬ 
lish  the  basic  formulas  to  be  used  below  and  in  succeeding  reports  that  will 
treat  acoustic  scattering  by  individual  localized  turbules  (or  vortons)  and  by 
ensembles  of  such  turbules. 

Section  2.1  contains  a  derivation  of  a  complete  set  of  basic  fluid  equations 
for  an  ideal  gas.  In  section  2.2,  we  assume  a  stationary  turbulence  model,  de¬ 
rive  equations  relating  the  turbulent  fields,  and  then  derive  the  desired  linear 
acoustic  differential  equations.  We  then  show'  that  Monin's  equation  results 
in  the  first-order  limit.  In  section  2.3,  we  develop  the  general  expressions  for 
the  Born  approximation  scattering  amplitudes  and  cross  sections  predicted 
by  the  equations  considered.  We  show  that  the  length  scale  of  the  scat¬ 
tering  volume  is  of  central  importance  in  predicting  the  general  behavior  of 
the  acoustic  scattering  at  wavelength  A  by  turbulence  of  scale  length  a.  For 
example,  we  show  that  the  short- wavelength  limit  of  Monin’s  equation  should 
be  used  only  if  both  a  >  A  and  a  >  a^,  and  that  the  standard  predictions 
from  Monin’s  equation  are  valid  for  any  ratio  a/A  only  if  a  <C  ag.  We  also 
obtain  the  interesting  general  result  that  the  first  Born  forward  and  back¬ 
ward  scattering  due  to  solenoidal  turbulent  flow  velocity  is  essentially  zero  if 
a  ttg,  but  nonzero  if  a  >  o^. 

In  contrast  with  the  turbulence  model  presented  here,  standard  treatments 
of  the  scattering  of  sound  weaves  by  atmospheric  turbulence  have  relied  on 
describing  the  turbulence  entirely  in  terms  of  the  autocorrelation  or  struc¬ 
ture  functions  of  the  turbulent  temperature  and  velocity  fields  [4,7,10].  This 
approach  has  been  particularly  successful  for  isotropic,  homogeneous,  fully 
developed  turbulence.  It  has  also  been  used  to  describe  scattering  by  inho¬ 
mogeneous  and/or  anisotropic  turbulence  [12-14].  Turbulence  has  also  been 
described  as  collections  of  individual  localized  eddies  or  vortices,  often  called 
turbules  [15,16].  Section  3  describes  research  aiming  at  the  description  of 
scattering  and  propagation  of  sound  waves  in  anisotropic  inhomogeneous  at¬ 
mospheric  turbulence. 


Section  3  considers  acoustic  scattering  by  individual  turbules.  Specifically,  we 
describe  a  general  scalable  localized  static  turbule  model,  valid  for  solenoidal 
turbule  velocities  and  arbitrary  morphology.  In  section  3.1,  we  use  the  results 
of  section  2  to  obtain  general  expressions  for  the  acoustic  scattering  ampli¬ 
tudes  and  cross  sections  due  to  scattering  by  the  velocity  and  temperature 
variation  fields  of  such  a  turbule,  in  first  Born  approximation.  Then  in  sec¬ 
tion  3.2,  we  describe  general  expressions  for  the  isotropic  ensemble  averages 
of  these  cross  sections.  Criteria  of  comparability  for  ensembles  having  differ¬ 
ent  envelope  functions  are  specified;  namely,  they  should  have  the  same  size 
and  the  same  kinetic  and  thermal  energy^  content.  We  present  two  examples, 
with  Gaussian  and  exponential  envelopes,  respectively,  and  contrast  their 
differential  and  total  scattering  cross  sections.  In  section  3.3,  we  also  derive 
the  azimuthally  dependent  velocity  scattering  cross  section  for  a  simple  ro¬ 
tating  turbule  with  a  spherically  symmetric  Gaussian  envelope.  In  addition, 
we  prove  a  general  theorem:  for  many  applications,  an  ensemble  of  randomly 
oriented  turbules  of  arbitrary  morphology  may  be  replaced  by  an  “equiva¬ 
lent”  ensemble  of  simple  turbules,  each  of  which  has  spherical  symmetry  but 
rotates  nonuniformly  about  a  fixed  but  randomly  directed  axis.  In  section  4, 
we  briefly  summarize  and  discuss  our  results. 
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2.  Basic  Theory 


2.1  Fluid  Equations 

In  descriptions  of  acoustic  propagation,  the  atmosphere  is  most  often  treated 
as  an  ideal  gas.  That  is,  viscosity  and  thermal  conductivity  are  neglected; 
the  fluid  is  assumed  to  be  in  local  thermal  equilibrium.  The  fundamental 
fluid  equations  for  such  a  gas  are  the  mass  conservation  equation,  the  Euler 
equation,  the  ideal  gas  equation  of  state,  and  the  thermal  energy  balance  or 
heat  flow  equation  without  the  heat  conduction  term: 


Ap  +  pV-v  =  0,  (1) 

Av +  -  g  =  0,  (2) 

p  =  pksT/M,  (3) 

DtU  +  (p  +  U)V  •  V  =  0.  (4) 


Here,  p(r,  t)  is  the  mass  density,  p(r,  t)  is  the  pressure,  r(r,  t)  is  the  temper¬ 
ature,  U{r,t)  is  the  internal  energy  density,  v(r,  t)  is  the  fluid  velocity,  /cg  is 
Boltzmann’s  constant,  g  is  the  vector  of  acceleration  of  gravity,  and  M  is  the 
average  molecular  mass,  which  is  constant  for  the  assumed  homogeneously 
mixed  atmosphere.  The  operator  A  is  the  convective  derivative, 

Dt  =  -|-  V  ■  V.  (5) 


The  heat  flow  equation  does  not  seem  to  have  been  used  in  previous  treat¬ 
ments  of  acoustic  propagation. 

For  a  fluid  in  local  thermal  equilibrium. 


AT 
2  ’ 


(6) 


where  F  is  the  number  of  degrees  of  freedom  per  molecule  excited  at  tem¬ 
perature  T.  We  have 


c,.  =  {dU/dT)p  =  pk^F/2M,  Cp  =  Cy  +  pk^/M,  7  =  Cp/cy  =  1  -h  (2/F). 


(7) 


Combining  equations  (3),  (4),  (6),  and  (7)  yields 

Dtp  +  7pV  •  V  =  0,  (8) 

Ar+(7-l)T'V-v  =  0,  (9) 

as  alternative  forms  of  the  heat  flow  equation.  Furthermore,  combining  equa¬ 
tions  (1)  and  (8)  yields 

Dt\n{p/ p^)  =  0,  (10) 

which  is  clearly  the  equation  for  isentropic  flow,  as  given  by  Tatarskii  [7]. 
Thus,  we  could  reverse  our  derivation  and  start  from  equations  (1-3)  and 
(10).  which  would  then  yield  equation  (8).  The  set  of  equations  (1-3)  and 
(8)  or  (1-3)  and  (10)  is  complete,  since  it  includes  six  equations  for  the  six 
flelds  {p.,p,  T,  Vi). 

2.2  Acoustic  Equations 

2.2.1  Linear  Acoustics 

For  linear  acoustics,  the  standard  approach  is  to  write  each  field  as  a  sum  of 
a  turbulent  part  and  an  acoustic  part,  and  retain  terms  in  the  resulting  equa¬ 
tions  only  to  first  order  in  the  presumably  small  acoustic  part.  Because  the 
time  scales  of  the  turbulent  flow  are  much  longer  than  those  of  the  acoustic 
fields,  the  time  dependence  of  the  former  is  often  neglected;  we  do  that  here. 
We  define 


=  Po(r)(l  +  c(r,  O); 

(11) 

=  Po{r){l+v{r..t)), 

(12) 

T(r,t) 

=  2^o(r)(n- 5(r,t)), 

(13) 

v(r,t) 

=  vo(r) +  u(r,t), 

(14) 

where  the  fields  with  subscripts  o  represent  the  turbulent  flow,  and  unsub- 
scripted  fields  represent  the  acoustic  flow.  We  also  take  the  turbulent  flow  to 
be  solenoidal, 

V  •  Vo  =  0,  (15) 

as  is  usually  done  [7].  Then  substitution  of  equations  (11-15)  into  equa¬ 
tions  (1-3)  and  (8)  yields  both  zero-order  equations  relating  the  turbulent 
fields,  and  acoustic  equations  linear  in  {e.p.S.u).  The  zero-order  equations 
are 

Vo  ■  Vpo  =  0,  Vo  •  Vpo  =  0,  po  =  Po^bTo/M,  (16) 
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(17) 


Vo  •  Vvo  +  Po  ^^^^0  -  g  =  0. 

The  acoustic  equations  are 

5t6  +  vo'V€  +  u-  Vlnpo  +  ^'U  =  0,  (18) 

diV  +  Vo  •  Vjy  +  u  •  V  Inpo  +  7^  ■  u  =  0,  (19) 

dtxi  +  Vo  •  Vu  +  u  •  Vvo  +  'y~^cl'Vrj  +  5Vpo/Po  =  0)  (20) 

6  =  ri-e,  (21) 

where  we  have  defined  the  adiabatic  sound  speed  Co(r)  by 

cl  =  7Po/po  =  (22) 

Combining  equations  (18)  and  (19)  yields  an  alternative  equation  relating 
rj.  e.  and  u: 

(at  +  vo-V)(7?-7e)  +  u- Vln(po/Po)  =  0.  (23) 


This  shows  clearly  that  ??  7^  7e  unless  po/po  is  uniform,  which  is  not  the  case 
in  general. 

The  acoustic  equations  (18—21)  are  a  complete  set  for  the  acoustic  fields 
(€,r?,  5,  u),  provided  (To(r),  vo(r))  are  given  and  equations  (16)  and  (17) 
and  boundary  conditions  are  used  to  obtain  (po(i'),  Po(r))- 

In  what  follows,  we  describe  situations  in  which  the  turbulence  is  localized 
'  in  some  bounded  volume  Vj--  Outside  Vr-  the  pressure,  temperature,  and 

mass  density  take  on  constant  uniform  values  (poo,  Tqc,  Poo),  the  flow  velocity 
Vo  =  0,  and  the  adiabatic  sound  speed  Coo  is  defined  by 

4  =  IPoo/Poo  =  'fhToo/M.  (24) 

2.2.2  First-Order  Acoustic  Equations 

In  the  atmosphere,  vq/coo  and  [To  -  Tod /Too  are  usually  of  the  same  order, 
both  much  less  than  unity,  so  it  is  a  good  approximation  to  discard  terms  from 
the  acoustic  equations  that  are  second  order  or  higher  in  these  ratios.  The 
principal  effect  of  gravity  in  the  atmosphere  is  to  produce  internal  gravity 
waves,  which  we  neglect.  Thus  equation  (17)  shows  that  Vpo  is  a  second- 
order  quantity.  Dropping  terms  containing  Vpo  from  fha  acoustic  equations 
(18-20)  leaves  only  two  relevant  coupled  equations,  namely,  the  first-order 
equations  (19)  and  (20).  Putting  in  time  dependence  exp(-  iut)  yields 

iuT]  =  7V-U  +  VO-V77,  (25) 

iujM  =  7“^CgVp -f- Vo  •  Vu -f  u  •  Vvq.  (26) 
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It  is  simple  to  eliminate  u  from  these  equations  and  get  a  wave  equation  for 
77  valid  to  first  order  in  uq/Cqc-  First,  replace  u  by  {iuj'f)  the  last 

two  terms  in  equation  (26) :  this  is  valid  to  first  order .  Then  evaluate  'y'V  •  u 
and  substitute  it  in  equation  (25).  This  procedure  yields 


+  =  V  •  [(Aro/roo)Vr7] 

+  (ifc/Coo){— Vo  •  V77  +  •  [vo  •  VV77  +  Vt?  •  Vvo]},  (27) 

where 

k  =  c’/coc,  ATo/Toc  =  1  -  To/Too  =  1  -  cl/cl^.  (28) 

Rearranging  terms  on  the  right  hand  side  of  equation  (27),  replacing  V^r]  by 
—  k'^rj  in  those  terms,  and  using  equation  (15)  yields  Monin's  equation  for  rj, 
written  here  in  summation  notation  with  di  —  d/dxi'. 

(V^  +  k‘^)7i  =  di[{ATo/T^)d,ri\  +  2iuj-^di{vojdjd^ii).  (29) 

This  is  the  accepted  wave  equation  for  the  acoustic  pressure  variation  that  is 
valid  to  first  order  in  the  small  stationary  quantities  {vq/c^,  ATo/Too).  The 
same  equation  is  valid  to  this  order  for  p'  =  pop,  since  Vpo  is  second  order. 
This  first-order  equation  is  valid  for  all  turbulent  length  scales  a.  If  these 
length  scales  are  all  very  large  compared  to  the  acoustic  wavelengths  [ka 
1).  then  derivatives  of  ATq  and  vq  may  be  neglected;  then  equation  (29) 
reduces  to 

(V^  +  k‘^)p  «  -  k^{ATo/Toc)p  -  {2ik/coc)vojdjP,  (30) 

which  is  the  same  as  the  first-order  wave  equation  given  by  Tatarskii  [4]. 

Note  that  our  formulation  uses  a  constant  reference  temperature  as  the 
temperature  outside  Vr-  Conventional  treatments  use  the  average  tempera¬ 
ture  T  in  Vr  as  the  reference,  which  might  or  might  not  be  equal  to  Toe- 


2.3  Scattering 

2.3.1  General  Formulation 

Equations  (29)  and  (30)  are  in  the  generic  form 

(V^  -)-  fc^)r?(r)  =  -  47r5(r)r7(r)  (31) 

where  ^(r)  is  an  operator  that  is  nonzero  only  inside  a  bounded  volume  Vr, 
or  goes  to  zero  rapidly  outside  this  volume.  The  operator  S  may  contain 
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derivatives  that  operate  on  This  standard  form  allows  a  causal  Green  s 
function  implicit  solution  for  Ti{r),  given  an  incident  plane  wave 

=  e'"',  (32) 

or  any  other  incident  wave  that  satisfies  the  Helmholtz  equation.  For  an 
incident  plane  wave,  the  implicit  solution  is 

r]{ri)  =  exp(ik  •  ti)  +  jd^r2ri2  exp(i/i:ri2)S'(r2)r?(r2),  (33) 

where  the  scattered  wave  r]s{Ti)  is  given  by  the  integral.  Here,  and  in  what 
follows,  ri2  =  ri  —  r2,  and  three-dimensional  spatial  integrals  with  no  limits 
indicated  involve  integration  over  all  space.  Actually,  the  integration  need 
extend  only  over  the  support  Vr  of  5(r2). 

In  the  far  field,  the  scattered  wave  has  the  form 

77,(r)  =  (34) 

where  the  scattering  amplitude  /(f)  for  the  scattering  volume  V'  is  given  by 

/(f)  =  f  d^riexp{-ikr  ■ri)S{ri)r]{ri),  (35) 

where  f  =  r/r  is  the  radially  outward  unit  vector. 

The  scattering  volume  V/  is  that  portion  of  the  volume  Vr  that  is  observed 
by  the  detector.  In  this  report,  we  assume  that  the  detector  is  in  the  far  field 
of  V'.  The  criteria  for  this  are 

rXOVA,  r»A,  (36) 

whichever  is  greater,  where  a'^  is  the  length  scale  of  Vg  and  A  =  27c/k  is 

the  acoustic  wavelength.  In  actual  experimental  situations,  these  criteria  are 
seldom  fulfilled. 

The  differential  and  total  scattering  cross  sections  ((T(f),  as)  of  the  turbu¬ 
lence  are  given  respectively  by 

<'{f)  =  l/(f)P.  =  Jdna(T)  (37) 

where  the  integration  extends  over  47r  sr  (solid  angle). 


2.3.2  Born  Approximation 

Ill  general,  the  correct  acoustic  wave  equation  would  involve  an  operator 

5  =  (38) 

where  Sn  is  order  in  the  small  quantities  {vofcoo,  ^Tq/Too)-  However, 
the  wave  equations  (29)  and  (30)  include  only  Si.  so  it  makes  no  sense  to  go 
beyond  the  first  Born  approximation  in  solving  for  the  scattering  amplitude 

/(r). 

Scattering  amplitudes.  We  get  the  first  Born  scattering  amplitude  by 
inserting  the  incident  (plane)  wave  (eq  (32))  in  the  integrand  of  equation  (35): 

/(f)  =  /dViexp(-i/:f  •  ri)5i(ri)exp(5’k-ri).  (39) 

J  Vs 

Here,  the  volume  Vg  is  now  determined  by  the  intersection  of  the  portion  V/ 
of  Vt  observed  by  the  detector,  and  the  portion  V7  of  V/  directly  illuminated 
by  the  source.  If  Vj  <  V/,  then  Vg  <  V'.  The  volume  K  might  be  regarded 
as  the  intersection  volume  of  two  cones,  the  detector  cone  of  view  and  the 
source  cone,  as  discussed  by  Tatarskii  [7].  For  convenience  in  what  follow's, 
we  ascribe  a  single  characteristic  length  scale  Ug  to  the  volume  Vg. 

For  wave  equation  (29),  the  operator  Si(r)  may  be  written 

Si{v)  =  SiT{r)  +  Si,{v),  (40) 

where 

S,T{r)  =  -j-^i({ATo/Too)^^^-{47^Too)-^[{^^ATo)^^  +  ATo^^^i],  (41) 

Att 

5i„(r)  =  -{i/27Tu)di{vojdjdi)  = -{i/2Truj)[{diVoj)djdi  +  vojdjdidi]  (42) 

are  the  operators  associated  with  the  turbulent  temperature  variation  and 
velocity,  respectively. 

It  is  clear  from  equation  (39)  that,  for  each  derivative  operator  dq  that  ap¬ 
pears  on  the  right  of  the  functions  ATq  or  vqj  in  ^i,  we  get  a  factor  ikg 
in  /(f).  However,  it  is  more  complicated  to  describe  the  general  effect  of  a 
derivative  operator  that  appears  on  the  left  of  these  functions.  It  must  be 
correct  to  substitute  equations  (41)  and  (42)  directly  into  equation  (39),  to 
get 

/(r)  = /T(f) +/v(r),  (43) 
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fj-{r)  =  {-k‘^/4TrToo)  j  oJ^ri  exp(- zK  •  ri)[z/c  •  ViAro(ri)  -  Aro(ri)],  (44) 

J  V^s 

=  (—  A:^/27rCoo)  f  exp(—  iK  •  ri)[—  ik~^k{kjdiiVQj{Ti)  +  fcjUoj(ri)],  (45) 

J  V5 

where 

K  =  kY-  k.  (46) 

But  then  we  cannot  draw  conclusions  about  the  general  behavior  of  / (f )  that 
are  independent  of  the  detailed  morphology  of  the  turbulence.  On  the  other 
hand,  it  must  also  be  correct  to  integrate  by  parts  in  equation  (39)  for  /(f). 
For  the  Sit  of  equation  (41),  this  yields 

frj.{Y)  =  {k‘^/4TrTc^)kih  [  AiAro(ri)  exp(- iK  •  ri) 

^  V5 

-{i/4nToo)ki  <{  dSi[AT’o(ri)  exp(- iK  •  ri)]^,,  (47) 

t/  Ss 

where  the  surface  integral  extends  over  the  surface  T>s  bounding  the  volume 

K. 

Consider  three  cases,  as  follows: 

Case  i.  All  the  length  scales  a  of  the  localized  turbulent  eddies  that  are  totally 
or  partially  included  in  Vs  are  much  less  than  dg-  This  is  usually  the 
case  for  inhomogeneities  in  the  inertial  range. 

Imagine  that  a  turbulent  eddy  or  turbule  is  localized  in  a  volume  “cen¬ 
tered”  about  a  random  location  b.  Then,  for  a  Os,  the  fields  (ATq,  uq)  of 
that  turbule  are  zero  on  the  surface  for  almost  all  b.  Thus,  the  surface 
integral  in  equation  (47)  is  negligibly  small  compared  to  the  volume  integral, 
and  the  scattering  amplitude  is  given  by  just  the  first  term  of  equation  (47) 
to  good  approximation: 

/r(f)  =  {k‘^/^T^Too)hri  [  dVi Aro(ri)  exp(-  iK  ■  ri).  (48) 

J  V5 

This  is  proportional  to  hfi  =  cosO,  where  6  is  the  scattering  angle.  This 
result  is  valid  for  all  ka,  for  a  <  a^.  The  corresponding  velocity  scattering 
amplitude  for  this  case  is 

/^(f)  =  -  {k'^/27rCoc){kiri)kj  f  dhiVoj{Ti)  exp(-  iK  •  ri),  (49) 

J  V5 
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which  also  displays  proportionality  to  cos^.  This  neglect  of  the  surface  inte¬ 
grals  is  usual  in  standard  treatments,  such  as  that  by  Tatarskii  [7].  Scattering 
amplitudes  and  cross-section  results  based  on  equations  (48)  and  (49)  have 
previously  been  given  [17]  for  a  particular  temperature  distribution  and  a 
particular  velocity  distribution. 

Case  ii.  All  the  length  scales  a  of  the  turbulent  eddies  in  Vr  are  greater  than 
as,  and  »  A,  so  /ca  >  1.  This  is  usually  the  case  for  large-scale 
inhomogeneities. 

In  this  case,  the  surface  integral  in  equation  (47)  is  not  negligible;  but  we 
can  neglect  {diATo,  diVoj)  in  equations  (41)  and  (42).  This  is  equivalent  to 
using  wave  equation  (30)  rather  than  equation  (29).  Dropping  the  first  term 
in  the  integrand  of  equation  (44),  we  get 

Mr)  ^  (kyAnToo)  /  dVi ATo(ri)  exp(-  fK  •  ri)  (50) 

J  V5 

to  good  approximation.  This  is  the  same  as  equation  (48),  except  for  a  factor 
of  cos  6. 

The  corresponding  velocity  scattering  amplitude  for  this  case  is  to  good  ap¬ 
proximation 

fv{r)  «  -  {k‘^/2nCoc)kj  [  dViUoj(ri)  exp(- iK  ■  ri),  (51) 

J  V^s 

which  is  the  same  as  equation  (49)  except  for  a  factor  of  cos  0. 

Case  iii.  A  >  a  >  a^,  so  A:a  <  1,  /ca^  <  1.  This  case  occurs  when  the  detector 
has  a  limited  field  of  view  and  turbule  sizes  are  small  because  of 
interference  of  the  ground. 

In  this  case,  the  surface  integral  in  equation  (47)  is  not  negligible  compared  to 
the  volume  integral,  and  neither  is  the  first  term  in  equation  (44)  compared 
to  the  second.  One  would  thus  expect  a  dip  in  the  scattering  near  9  =  90°,  but 
not  a  zero;  the  scattering  amplitude  /r(r)  of  equation  (44)  directly  involves 
the  derivative  of  ATq  (which  it  does  not  for  cases  (i)  and  (ii)  above)  and 
similarly  for  The  full  expressions  (44)  and  (45)  should  be  used. 

It  seems  worthwhile  to  note  that,  if  a  ':$>  ag,  then  both  ATo(r)  and  vo(r)  do 
not  change  appreciably  within  K;  a  turbulent  eddy  with  such  a  large  scale 
length  appears  to  the  observer  to  be  a  uniform  wind  field  and  a  uniform 
temperature  shift  over  the  volume  Vg.  For  such  cases,  the  fields  Aro(ri)  and 
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Voj  may  be  taken  outside  the  integrals  in  equations  (44),  (45),  (50),  and 
(51),  and  evaluated  at  any  point  inside  14,  say  at  some  center  point  r*.  For 
example,  for  very  large-scale  inhomogeneities,  with  a  dg,  ka  ^  1,  equation 
(50)  reduces  to 

f^(r)  =  {kyA7rT^)ATo{Ts)  [  dVi  exp(-  iK  •  ri).  (52) 

In  summary,  the  scale  length  of  the  scattering  volume  is  at  least  as  impor¬ 
tant  as  the  scale  lengths  a  of  the  turbulent  eddies  and  the  acoustic  wavelength 
A  =  27r/iv  in  determining  the  general  behavior  of  the  first  Born  acoustic  scat¬ 
tering  amplitudes.  If  a  <  a*,  then  for  any  ka  values,  equations  (48)  and  (49) 
result.  If  a  >  a,  and  ka  »  1,  then  equations  (50)  and  (51)  result,  which 
are  the  same  as  equations  (48)  and  (49),  except  for  the  absence  of  the  cosd 
factor,  and  the  same  as  would  be  achieved  directly  if  we  started  from  the 
approximate  wave  equation  (30)  instead  of  from  equation  (29).  If  o  >  as, 
but  the  ka  are  not  large,  then  the  complete  equations  (44)  and  (45)  must  be 
used. 

Cross  sections.  From  equations  (37)  and  (43),  the  first  Born  differential 
cross  sections  are  given  by 

cr(f)  =  (Jr(f)  +  <7’t;(f)  -f-  crTv(f’)  (53) 

where 

Mr)  =  |/T(f)P,  (Tvir)  =  |/t;(f)|^  aTv{r)  =  /r(r)/;(f) -i- c.  c.  (54) 

If  we  are  interested  in  the  ensemble  average  over  assumed  stochastic  ho¬ 
mogeneous  turbulent  fields  (ro(r),  Vo(r)),  and  we  assume  that  fluctuations 
in  ro(r)  are  uncorrelated  to  fluctuations  in  vo(r),  and  (vo(r))  =  0,  then 
((^Tt;(f))  =  0. 

Consider,  for  example,  the  case  for  (o-T(f))  for  a  <  a*.  From  equations  (54) 


and  (48), 

K(f))  =  (J:Vl67r^r|,)cos^«{|Afo(K)p) 

(55) 

where  , 

Afo(K)  =  /  dVi Aro(ri)  exp(-  iK  ■  ri) 

JVs 

(56) 

is  the  Fourier  transform  of  ATo(ri)  with  respect  to  the  volume  14 .  Then  for 
homogeneous  turbulence, 


(|Aro(K)p)  =  /  d^ri  f  dV2Br(ri2)exp(-iK  ■  ri2),  (57) 

J  Vs  Vs 


where 


Brirn)  =  {ATo{r,)ATo{v,))  (58) 

is  the  autocorrelation  function  of  the  turbulent  temperature  variations.  Since 
a  <§:  as,  the  scale  length  of  Brirn)  is  also  much  less  than  a*.  That  is,  Brirn) 
goes  to  zero  for  ri2  >  some  a',  where  a'  Oj.  Therefore,  the  integration 
variables  (ri,  r2)  in  equation  (57)  may  be  changed  to  r  =  ri2,  R  =  (ri+r2)/2, 
and  the  limits  on  the  integral  over  r  may  be  changed  to  all  space,  while  the 
integral  over  R  remains  over  14,  with  negligible  error.  Thus  we  get 

(ar(f ))  =  {kyieTT^Tl)  cos2  eVsBriK)  (59) 

where 

.^^(K)  =  J  d^rBr(r)exp(-iK  •  r)  (60) 

is  the  full  Fourier  transform  of  the  autocorrelation  ^^(r)  integrated  over  all 
space,  to  good  approximation.  This  formula  is  generally  used  in  conventional 
treatments  of  homogeneous  turbulence  [7]. 

However,  consider  the  opposite  limits,  a  :$>  Ug,  ka  ^  1.  Then  from  equa¬ 
tions  (54)  and  (52), 

(ar(r)}  {k'^ /16-k‘^T^)Bt{0)\  (fn  exp(-  iK  ■  ri)p.  (61) 

This  is  not  only  missing  the  cos^  6  factor,  but  it  also  has  a  quite  different 
dependence  on  K  from  that  in  the  expression  in  equation  (59). 

In  the  companion  report  [2],  we  will  assume  that  14  is  large  enough  so  that 
a  Qs  for  all  turbulent  scale  lengths  of  interest,  so  that  equations  (48)  and 
(49)  apply.  This  will  allow  us  to  compare  results  using  ensembles  of  local¬ 
ized  turbulent  eddies  or  turbules  of  different  scale  lengths  with  the  standard 
structure  function  results. 

Forward  and  backward  velocity  scattering.  Consider  equation  (49)  for 
the  Born  scattering  amplitude  due  to  turbulent  velocity  in  case  (i),  a  <C  o*. 
In  the  forward  direction  f  =  k  or  K  =  0,  we  have 

fyCk)  =  {k'^/2TTCoc){cos6)kj  f  dVuoj(r).  (62) 

J  Vs 

Consider  the  following  integral: 

Ij  =  (PrdiivoiTj)  =  d^rvoj{r).  (63) 
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Since  the  velocity  is  solenoidal,  diVQi  =  0.  But  by  the  divergence  theorem, 

/,•  =  f  dEiivoiVjh,.  (64) 

By  the  same  reasoning  used  earlier  for  a  a^,  this  surface  integral  vanishes 
(is  negligibly  small)  because  for  almost  all  positions  of  the  center  of  any 
localized  eddy,  voi  is  zero  on  the  bounding  surface  Therefore,  if  all  a  < 
as-.  ^  0. 

Consider  the  Fourier  transform  of  Uoi(r)  in  the  volume  Vg'. 

voi{K)  =  f  d^re~^^'''voi{r).  (65) 

J  V5 

Again,  because  diVQi  =  0,  we  have 

KiVoi{K)  =  if  d^rdi[voi{v)e~"^''']  -  i  f  dEi[voi(i')e"''**']E,-  (66) 

JVs 

For  the  same  reason  as  before,  the  surface  integral  is  vanishingly  small  if 
a  C  Os-  From  equation  (49),  the  Born  scattering  amplitude  in  the  backward 
direction  r  =  -  k  is  given  by 

/„(-k)  =  -{k'^/27rCoo){cosB)kiVoi{-2k),  (67) 

and  is  therefore  zero  (negligibly  small).  This  zero  result  for  backward  velocity 
scattering  was  noted  by  Tatarskii  [7]. 

On  the  other  hand,  for  cases  (ii)  and  (iii),  a  >  a*,  the  surface  integrals 
(64)  and  (66)  are  not  negligibly  small,  and  the  results  /v(±k)  «  0  do 
not  hold.  In  fact,  for  a  »  Ug,  we  have  kjlj  VgkjVoj{rg),  and  KiVoi{K)  = 
*^oi(rs)  /s'  Both  are  nonzero  in  general.  This  is  what  would  be 

expected  if  we  were  observing  a  volume  Vg  in  which  there  is  a  uniform  wind 
field. 


3.  Scattering  from  Model  Turbules 


In  section  2,  we  describe  some  general  properties  of  acoustical  scattering  by 
atmospheric  inhomogenieties,  which  are  true  regardless  of  the  internal  con¬ 
figuration  or  morphology  of  the  inhomogeneity.  In  this  section,  we  describe 
scattering  properties  of  inhomogeneities  with  prescribed  morphology,  as  well 
as  properties  of  isotropic  ensembles  of  inhomogeneities  with  similiar  mor¬ 
phology. 

We  define  a  turbule  as  a  localized  eddy  or  vortex,  characterized  by  a  location 
b,  a  scale  length  a,  and  flow  velocity  and  temperature  fields  (vq,  ATq).  Here, 
“localized”  means  that  these  fields  go  to  zero  rapidly  for  |r  -  b|  >  a. 

In  actuality,  the  fields,  the  location,  and  the  scale  length  depend  on  time. 
As  is  often  done  in  discussions  of  interactions  of  acoustic  waves  with  turbu¬ 
lence,  we  neglect  this  time  dependence,  and  instead  treat  the  location  and 
the  morphological  parameters  in  (vq,  ATq)  as  stationary  but  stochastic  vari¬ 
ables.*  The  fields  (vq,  ATo)  are  assumed  independent,  and  may  have  arbitrary 
morpholog>^ 

We  assume,  as  usual,  that  the  turbule  flow  is  solenoidal: 

V-vo  =  0  — >  Vo  =  Vx  Aq.  (68) 

Here,  Aq  is  a  vector  field. 

Since  we  are  interested  ultimately  in  ensembles  of  self-similar  turbules  of 
different  scale  lengths,  we  assume  that  the  functional  forms  of  (Aq,  ATq) 
satisfy 

Ao(r)  =  aA(^),  Aro(r)  =  Too  -  To(r)  =  T(^),  where  ^  =  (r  -  b)/a.  (69) 

Here,  To(r)  is  the  temperature  field,  and  Too  is  the  uniform  constant  reference 
temperature  outside  the  turbule.  These  are  general  scalable  forms,  chosen  for 
convenience  so  that  the  functions  (A(^),T(^))  have  the  dimensions  (velocity, 
temperature),  respectively.  Then  equation  (68)  yields 

vo(r)  =  X  A(^).  (70) 

*This  means  that  Doppler  shifts  cannot  be  described. 
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3.1  General  Expressions 

We  use  the  expressions  for  the  first  Born  scattering  amplitudes  derived  in 
section  2,  equations  (48)  and  (49),  for  a  <C  a^,  where  is  the  scale  length  of 
the  scattering  volume: 

/^(r)  =  -(fcV47rr^)cos^Afo(K),  (71) 

fv{r)  =  -  (A:^/27rCoo)cos0  k  •  vo(K).  (72) 

Here,  (Afo(K),  vo(K))  are  the  Fourier  transforms  of  (Aro(r),  vo(r)),  respec¬ 
tively,  and  for  acoustic  wave  time  dependence  exp{—iu)t), 

k  =  io/coo,  K  =  k{r  -  k),  K  =  2k  sin(6*/2),  (73) 

where  6  is  the  (polar)  scattering  angle,  that  is,  the  angle  between  the  incident 
plane-wave  propagation  direction  k  and  the  observation  direction  f .  Also, 

Coo  =  (7^-bToc/M)5  (74) 

is  the  adiabatic  sound  speed  outside  the  turbule;  here,  7  is  the  specific  heat 
ratio,  M  is  the  mean  molecular  mass,  and  k^  is  Boltzmann's  constant.  The 
Fourier  transforms  are  defined  by 

Afo(K)  =  ydViAro(ri)e'**^'’b  xo{K)  =  (75) 

where  the  integrals  here  have  been  extended  over  all  space,  since  a  <C  ag. 

Combining  equations  (68-75)  yields  the  following  general  expressions  for  the 
scattering  amplitudes  for  a  turbule  localized  around  r  =  b: 

/r,.(f,b)  =  e-*77-..(f).  (76) 

where 

f^(r)  =  -(itV/47rroo)cos0f(Ka),  (77) 

f^{r)  =  -i(/c^a^/27rcoo)sin0cos0[<^  •  A(Ka)]  (78) 

are  the  scattering  amplitudes  for  a  turbule  localized  around  the  origin. 

The  scaled  Fourier  transforms  (f ,  A)  are  given  by 

f(y)  =  A{y)  =  (79) 
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and 


=  —  ei  sin  ;^  +  02  cos  v?  (80) 

is  a  unit  vector  in  the  direction  of  increasing  azimuthal  scattering  angle  9?. 
The  Bi  (i  =  1,  2,  3)  are  Cartesian  basis  vectors. 

Note  the  sin  0  cos  ^  dependence  in  fy'.  this  occurs  for  arbitrary  turbule  mor¬ 
phology,  and  implies  that  the  first  Born  scattering  is ^ero  in  the  forward  and 
backward  directions,  as  well  as  at  6  =  7r/2,  smce  A(y)  is  bounded  every¬ 
where  for  a  localized  bounded  A(^).  That  fv{k)  —  /„(— k)  =  0  was  shown 
in  section  2  to  be  a  general  consequence  for  localized  solenoidal  turbulent 
velocities  for  a  <C  equations  (68-78)  constitute  another  general  proof  of 
that  result. 

Another  useful  expression  for  fy  results  from  using  equation  (73)  in  equa¬ 


tion  (78): 

fvir)  =  —  i{k'^a^ /2'iTCao)  cos 6 cos{6f  2) {K a) {ip  ■  A(Ka)).  (81) 

The  differential  scattering  cross  sections  are  given,  as  usual,  by 

cr(f)  =  |/(f,b)p,  /(f,b)  = /v(f,b) +  /r(r,b).  (82) 

Thus 

a(r)  =  crr(r)  -I-  (7„(f)  -h  (Ttv{^)  (83) 

where 

ar{T)=\Mr)\^  =  (/cV/lGTr^T^)  cos^  (Ka)|2  ^  (84) 

^.(f)  =  l/v(f)l"  =  {k^ay47r^cl)cos^dcosHe/2){Ka)^\>p-A{Ka)\\  (85) 

(TTvir)  =  /T(f)/;(f)-l-c.  c.  (86) 


Note  that  these  cross  sections  are  independent  of  turbule  location  b.  Also 
note  the  implicit  dependence  of  ay  on  the  azimuthal  scattering  angle  p. 

The  total  scattering  cross  sections  are  the  integrals  of  these  over  47r  sr  (solid 
angle). 

We  cannot  go  further  unless  we  choose  specific  functional  forms  for  (T(^), 
A(^)),  or  unless  we  average  general  forms  over  orientations  and  morphologies. 
In  section  3.2,  we  consider  isotropic  ensembles  of  turbules  of  scale  length  a.  In 
section  3.3,  we  choose  instead  one  simple  specific  form,  in  order  to  illustrate 
the  azimuthal  dependence  of  the  velocity  scattering  that  must  occur  (but 
which  disappears  for  isotropic  ensembles).  We  also  use  this  simple  form  to 
show  that  general  isotropic  ensembles  may  often  be  replaced  by  ensembles  of 
spherically  symmetric,  nonuniformly  rotating  turbules. 
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3.2  Isotropic  Ensembles 

In  a  follow-up  report  [2]  we  will  compare  (1)  the  scattering  produced  by 
isotropic  ensembles  of  turbules  of  different  scale  lengths  and  locations  with  (2) 
the  scattering  predicted  by  the  conventional  structure  functions  of  isotropic 
homogeneous  fully  developed  turbulence.  We  will  need  ensemble  averages 
over  the  orientation  parameters  of  the  individual  turbules;  thus  we  derive 
these  here. 

3.2.1  Orientation  Averages 

Imagine  a  static  turbule  of  scale  length  a,  such  that  in  a  “primed "  frame  of 
reference  (its  “rest  frame”),  it  has  arbitrary  morpholog}^.  A  subensemble  is 
formed  by  the  application  of  different  static  rotations  R{ip)  to  this  generic 
form,  with  all  rotations  equally  probable  for  an  isotropic  case.  Here  xp  = 
(^i)'02,'03)  are  the  three  Euler  angles  that  characterize  a  rotation.  Consider 
first  the  temperature  field.  In  the  primed  frame,  we  may  write  in  general 

|r(y')P  =  EemaeMYimiO',  ^')  =  |r(y)P,  (87) 

where  the  last  equality  follows  because  of  scalarity,  and  the  sum  is  a  generic 
expansion  in  spherical  harmonics  in  the  primed  frame,  with  given  but  arbi¬ 
trary  aem{y)-  But  we  know  that  the  spherical  harmonics  in  the  two  frames 
are  related  by 

Ye^ie',  ^')  =  ^),  (88) 

where  the  Z)^^/('0)  are  the  irreducible  representations  of  i?(^)  for  integer 
For  all  rotations  equally  probable,  we  then  want 

=  (87r^)'^  d'lpz  =  (89) 

Thus, 

{|f  {y)p)^  =  a<x{y)  =  (90) 

which  is  a  function  of  y  =  |y|  only.  Here,  6T  is  a  temperature  variation 
amplitude  parameter.  We  do  not  specify  the  dimensionless  real-valued  enve¬ 
lope  function  Bxiy)  at  this  time;  we  present  two  examples  (Gaussian  and 
exponential)  in  section  3.2.3. 

Consider  now  the  velocity  field.  For  an  isotropic  ensemble,  each  component  of 
an  otherwise  unrestricted  vector  field  must  be  statistically  independent  and 
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have  zero  mean,  and  each  component  must  obey  the  same  statistics.  Thus 
we  may  put 

{My)h  =  0’ 

(ii(y)i*(y))^  =  (91) 

where  v  is  a  velocity  amplitude  parameter.  The  envelope  function  By{y)  may 
be  different  from  Briy)-  It  is  important  to  note  that  equation  (91)  yields  the 
usual  form  for  (f;oi(K)f^  *oi  (K))  =  ^ij(K)i  namely, 

$ij(K)  =  {S,j  -  KiKj)  (function  of  IK]).  (92) 

This  follows  because  of  equation  (70). 

From  equations  (84)  to  (86),  w^e  have,  for  an  isotropic  ensemble  of  scale  length 
a, 

^7’t;(f)  =  0,  (93) 

ar{r)  =  i7ra^){5T/4T^)HkaY  cos'^  0B^{Ka),  ^  (94) 

a,,(f)  =  \{T:a^){vl2c^f{ka)Uos^ecos^{ei2){KafBl{Ka).  (95) 

The  implicit  !/;-dependence  of  ay(r)  in  equation  (85)  for  an  individual  turbule 
is  eliminated  by  the  isotropic  averaging  process. 

3.2.2  Comparable  Turbules 

It  is  important  to  be  able  to  compare  the  scattering  produced  by  ensembles 
of  the  same  scale  length,  but  different  morphologies  (i.e.,  ensembles  char¬ 
acterized  by  different  envelope  functions  B(y)).  In  this  report,  we  choose 
comparable  turbules  to  be  those  having  the  same  rms  radius  and  the  same 
average  energy  content. 

One  way  to  define  an  rms  radius  is  to  write 


=  JdHeBii)/  fdHBd), 

(96) 

where  B{^)  is  the  inverse  Fourier  transform  of  B{y), 

B{i)  =  (21^)-^  jdSe‘^yB(v). 

(97) 

This  is  equivalent  to 

(<*,„/a)  =  [  -  VjS(y)/S(!,)]J,„. 

(98) 
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We  adopt  this  definition. 

The  kinetic  energy  content  of  the  ensemble  average  turbule  of  scale  length  a, 
to  lowest  order  in  Vq  (second  order),  is  proportional  to 

jd^ivUr)).  (99) 

Substitution  of  equations  (69),  (79),  and  (91)  shows,  after  a  little  algebra, 
that  this  integral  is  proportional  to 

=  i»=(27r)-V2  jdWSUv)-  (100) 

We  take  the  excess  internal  (thermal)  energy  content  of  the  ensemble  average 
temperature  turbule  to  be  proportional  to 

=  [(8/(2,r)=/“)  =  [(8/(2^)V^)  jd‘y{\f(y)\-‘}]K  (101) 


3.2.3  Examples 

(a)  Gaussian  envelope.  For  convenience,  we  let 

BT(y)  =  BM  =  B,[y)  =  e-«'l\  (102) 

Such  forms  have  been  used  in  the  literature  [15].  Equations  (98),  (101),  and 
(102)  then  yield 

flrmi/a  =  \/3/2,  Exg  ~  STg,  Eyg  =  v^.  (103) 

The  ensemble  average  differential  scattering  efficiencies  follow  from  equations 
(94)  and  (95): 

^rs(r)/7ra^  =  (STg / {ka)"^  cos^  6 exp  [- {ka -  cos 9)],  (104) 

at,g(f)/7ra^  =  i(up/2coo)^(/i;a)®sin^6>cos^6>exp  [  -  (A:a)^(l  -  cos6>)],  (105) 

where  we  define  a  scattering  efficiency  as  the  ratio  of  a  cross  section  to  tto^. 
These  differential  scattering  efficiencies  are  plotted  in  figures  1  and  2  as 
functions  of  scattering  angle  9  for  several  different  size  parameters  {ka),  nor¬ 
malized  to  unit  values  of  6Tg/Too  and  Vg/coo-  These  expressions  yield  exactly 
the  same  dependence  of  the  scattering  on  9,  k,  and  a  as  the  expressions  in 
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Figure  1.  Temperature  scattering  efficiency  versus  scattering  angle  9,  for  6T/Tx  = 

1.0,  for  a  Gaussian  isotropic  ensemble  of  turbules  of  scale  length  a,  for  several  size  param¬ 
eters  ka  -  27ra/A,  A  =  acoustic  wavelength,  in  Born  approximation. 


Scattering  angle,  0  (°) 

Figure  2.  Velocity  scattering  efficiency  a^{d)lT:o?  versus  scattering  angle  9,  for  V/Coc  =  1-0, 
for  a  Gaussian  isotropic  ensemble  of  turbules  of  scale  length  a,  for  several  size  parameters 
ka  =  2-KalX,  A  =  acoustic  wavelength,  in  Born  approximation. 
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Mellert  [18]  for  a  Gaussian  spectrum  of  temperature  and  wind  velocity  fluc¬ 
tuations.  That  is,  a  Gaussian  envelope  function  for  an  isotropic  ensemble 
of  turbules  of  a  single  scale  length  is  the  same  as  a  Gaussian  spectrum  for 
isotropic  turbulence  of  a  single  scale  length,  as  would  be  expected. 

The  total  efficiencies  are  given  by 

=  2Tv{6TglAToo?{x  -  2 +  2/x)  -  +  2  +  2/x).,  (106) 

=  (27r/3)(vp/2C(x,)^a:[(l -h  e“^'^)(2 -I- 24/x^) 

-  (1  -  e-2^)(10/a:  +  24/x^)],  (107) 

where  x  =  (fca)^. 

In  figure  3,  the  numerical  integration  of  the  differential  scattering  efficiencies 
are  plotted  versus  ka  for  unit  6T/T00  and  v/c^o-  For  small  ka.  the  total 
efficiencies  reduce  to 


GTg/'Xa^  ~ 

{^'xl'i){bTgl^T^)\ka)\ 

(108) 

a^gj-XQ^ 

\{^nl\^){vgl2c^f{kaf. 

(109) 

For  large  ka,  they  reduce  to 

(J  'I'g  f  TT  a 

2'K{6Tg/AToof{kaf.. 

(110) 

a  xjg  j  'xa 

{ATc/2,){vg/2cocf(kaf. 

(111) 

Figure  3.  Total  Born  approximation  (temperature,  velocity)  scattering  efficiency 
{Wt I'kc? versus  size  parameter  ka  =  27ra/A,  A  =  acoustic  wavelength,  a  =  turbule 
scale  length,  for  a  Gaussian  isotropic  ensemble,  for  6T/Toc  =  v/Coc  =  1-0. 


As  is  well  known,  the  Born  approximation  fails  in  the  limit  of  arbitrarily 
large  ka\  the  actual  total  scattering  efficiencies  should  be  independent  of  ka 
in  this  limit.  But  we  cannot  simply  apply  higher  order  approximations  for  the 
scattering  amplitudes  for  large  ka.  because  the  wave  equations  used  are  valid 
only  to  first  order  in  uo/coc-  In  order  to  obtain  trustworthy  expressions  for 
scattering  cross  sections  for  very  large  size  parameters,  we  must  first  obtain 
and  solve  wave  equations  valid  to  all  orders  of  Vq/coo  and  its  derivatives. 

Several  features  of  the  above  results  are  remarkable.  If  (ST/Too)  ~ 
then  Oyg  axg  for  small  ka,  while  a^g  ^  oxg  for  large  ka.  Also,  — 

0  for  ^  =  (0, 7r/2,  tt),  while  is  a  maximum  in  the  forward  direction. 

(b)  Exponential  envelope.  Again  for  convenience,  we  let 

Briy)  =  BM  =  B.(y)  =  (1  +  (112) 

where  o:  is  a  parameter  to  be  adjusted.  This  corresponds  to  the  following 
envelope  function  in  position  space: 

=  {2-k)-^  (x  (1  +  aOe-"^.  (113) 

The  presence  of  the  factor  (1  +  is  necessary  in  order  to  ensure  that  the 
turbule  velocity  field  is  bounded  at  ^  =  0. 

Equations  (98),  (100),  and  (101)  then  yield 

arms/a  =  \/l8/a,  Exe  = 

Comparing  equation  (103)  yields 

a  =  VT2,  6T^  =  (0.84)((5Tp),  =  0.64^^;  (115) 

we  make  this  comparison  so  that  these  exponential  model  turbules  of  scale 

length  a  have  the  same  ensemble  average  values  as  the  Gaussian  model  for 
{arms-,  Exe-.  l-e.,  SO  that  the  two  ensembles  of  different  average  morphology 
be  “comparable.” 

The  ensemble  average  differential  scattering  efficiencies  then  follow  from 
equations  (94)  and  (95): 

axe{r)lTta^  =  {0. 70) {6Te / cos'^  ^[l  +  I {ka)%l  -  cos  e)]~\  (116) 

ave{r)/TTa‘^  =  |(0.41)(ue/2coo)^(A':a)®sin^^cos^^[l  +  i(A:a)^(l  -  cos^)]  .  (117) 
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These  are  plotted  in  figures  4  and  5  as  functions  of  scattering  angle  6  for 
several  size  parameters  ka,  for  unit  values  of  STe/T^o  and  Ve/coo- 

The  calculations  for  the  total  cross  sections  versus  ka  are  algebraically  quite 
complicated.  We  give  the  small  ka  limits  only;  the  curves  of  figure  6  were 
found  by  direct  numerical  integration  of  the  differential  cross  sections,  equa¬ 
tions  (116)  and  (117).  For  small  size  parameters,  we  have  for  the  total  effi¬ 
ciencies 

are/Tta^  ~  (47r/3)(0.70)(<5r^/4roo)"(A:a)'  (0.70) (ar^/™'),  (US) 

~  5(0.41)(87r/15)(up/2coo)^(^fl)^  ~  (0.41)(crt,g/7ro  ).  (HO) 

Thus,  for  small  ka,  the  total  cross  sections  for  the  exponential  ensemble  are 
smaller,  but  of  the  same  order  as  those  for  the  Gaussian  ensemble.  Figures  3 
and  6  allow  a  comparison  of  the  total  efficiencies  versus  ka. 


Figure  4.  Temperature  scattering  efficiency  crTOIna^  versus  scattering  angle  9,  for 
ST/Toc  =  1-0,  for  an  exponential  isotropic  ensemble  of  turbules  of  scale  length  a,  for 
several  size  parameters  ka  =  27ra/A,  A  =  acoustic  wavelength,  in  Born  approximation. 


Figure  5.  Velocity  scattering  efficiency  a^eiiTO?  versus  scattering  angle  6.  for  v/coc  - 
1.0,  for  an  exponential  isotropic  ensemble  of  turbules  of  scale  length  a,  for  several  size 
parameters  ka  =  2Tia/X,  A  =  acoustic  wavelength,  in  Born  approximation. 


Figure  6.  Total  Born  approximation  (temperature,  velocity)  scattering  efficiency 
(uT/7ro2,a„/7ra2)  versus  size  parameter  ka  =  27ra/A,  A  =  acoustic  wavelength,  a  =  turbule 
scale  length,  for  an  exponential  isotropic  ensemble,  for  6T/Toc  =v/Coc  =  1-0. 
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3.3  Individual  Turbule  Model 

This  section  illustrates  the  azimuthal  dependence  of  the  velocity  scattering 
from  even  the  simplest  model  turbule,  and  shows  that  all  isotropic  ensembles 
may  indeed  be  composed  of  such  simple  model  turbules.  For  example,  suppose 
we  choose 

A(^)  =  ue“^^  (120) 

where  u  is  a  constant  vector,  as  the  simplest  Gaussian  turbule  model.  Then 
we  have,  from  equation  (69), 

vo(r)  =  X  A(^)  =  2u  X  (121) 


If  we  define 


u  =  iafi 


(122) 


and  use  the  definition  of  ^  (eq  (69)),  we  get 

vo(r)  =  [O  X  (r  -  b)]  exp  -(r  -  b)^/a^ 

Thus,  this  model  of  a  single  turbule  represents  a  nonuniform  rotation,  with 
angular  velocity  parameter  fl,  of  a  spherically  symmetric  Gaussian  turbule 


(123) 


of  scale  length  a,  centered  at  r  =  b. 

The  Fourier  transform  of  A(^)  is  given  by  equation  (79);  we  have 

A(y)  =  (124) 

We  write  u  in  terms  of  its  polar  and  azimuthal  angles: 

u  =  ei  sin  cos  ipu  +  ^2  sin  9u  sin  ipu  +  ^3  cos  6u-.  (125) 


whereby  from  equations  (80)  and  (85),  the  first  Born  differential  scattering 
efficiency  is 


ay{r)/7ia‘^  =  {ka)^{v/2c^f  sin^  Ocos^  6sm^  [-(^a)^(l  -  cos6')]  .  (126) 


This  exhibits  strong  azimuthal  dependence.  It  is  zero  if  =  0  or  7r(u  =  ±k), 
and  if(p  =  {(fui  <^ti±7r).  That  is,  it  goes  to  zero  if  u  is  in  the  plane  of  incidence, 
the  plane  whose  normal  is  k  x  f . 

If  we  average  over  random  directions  of  u,  using 

(47r)"^  /  difu  f  d^iisin0„[sin^6'usin^((/?u  -  V?)]  =  1/3,  (127) 

Jo  Jo 
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we  get 


at,(f)/7ra^  =  ^{kaf  {u/2coof  siri^  9 cos^  6 exp  [-{kaf{l  -  cos6')  ,  (128) 

which  is  identical  with  equation  (105)  for  the  Gaussian  ensemble,  if  u  =  Vg. 
In  fact,  from  equations  (124)  and  (127),  we  have 

(i.(y)-4'(y)>  =  (129) 

which  yields,  from  equations  (91)  and  (102),  for  the  choice  u  =  Vg, 

B2(y)  =  e-*''/^  (130) 

the  same  envelope  function  as  that  chosen  for  the  Gaussian  ensemble.  The 
same  conclusions  could  be  drawn  for  all  individual  turbules  of  the  generic 
form  _ 

A(0  =  u/(0  A(y)  =  u/(y),  (131) 

where  u  is  a  constant  vector  that  may  have  any  direction  and  magnitude, 
and  /  is  a  scalar  function  of  ^  =  |^|.  The  corresponding  velocity  field  is 

vo(r)  =  X  A(^)  =  -  (u  X  O(rV'(0)-  (132) 

If  we  define  ft  as  in  equation  (122),  and  /i(0  =  -  K"V'(0:  then  we  have 

vo(r)  =  X  (r  -  b)]/i(|r  -  b|/a).  (133) 

This  shows  that  the  vector  u  is  along  the  spin  axis  ft  of  any  simple  model 
turbule.  All  such  “simple"  turbules  are  rotating,  in  general  nonuniformly, 
with  spherically  symmetric  envelopes. 

For  random  orientations  of  u,  we  have  {uiUj)  =  and  thus 

(ii(y)i*(y))  =  iM^|/(y)P%,  (134) 

whereby  from  equation  (91),  we  can  identify 

BHy)  =  (135) 

as  the  corresponding  isotropic  ensemble  envelope  function. 

Similar  considerations  apply  for  the  temperature.  Consider  a  spherically  sym¬ 
metric  turbule  with  A7o(r)  =  T{^).X  =  |CI-  Then, 
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\ny)?  =  (If  (»)!")  =  ASTfB^(vh  (136) 

since  orientation  averaging  has  no  effect  on  a  spherically  symmetric  function. 

Based  on  these  considerations,  we  can  reach  the  following  important  con¬ 
clusion:  Every  ensemble  of  randomly  oriented  static  turbules  of  a  given 
scale  length,  within  which  each  member  has  arbitrary  morphology,  is  equiv¬ 
alent  to  a  single  turbule  with  a  spherically  symmetric  temperature  variation 
Aro(r)  =  T{^),  and  to  a  collection  of  simple  turbules  with  velocity  vector  po¬ 
tential  A(£)  (given  by  eq  (131),  with  u  a  randomly  oriented  constant  vector). 
The  specific  equivalence  is  established  by  equations  (135)  and  (136).  Here, 
“equivalent”  means  “produces  the  same  (incoherent)  acoustic  scattering.” 
This  theorem  is  similar  to  that  for  “equivalent  spheres”  that  may  replace  a 
collection  of  randomly  oriented  dielectric  particles  of  arbitrary  morphology 
in  incoherent  electromagnetic  scattering. 

As  in  the  electromagnetic  case,  the  theorem  is  valid  only  for  incoherent  scat¬ 
tering  from  an  actual  collection  of  turbules.  But  this  case  is  of  great  im¬ 
portance  and  generality.  On  the  other  hand,  if  coherence  is  important  in 
the  calculation  of  multiple  scattering  from  an  actual  collection  of  randomly 
oriented  turbules  of  general  morphology,  then  we  should  not  perform  the  ori¬ 
entation  averaging  until  after  we  have  solved  for  the  scattering  cross  sections 
of  the  whole  collection,  for  each  of  many  realizations  of  the  collection.  In 
such  cases,  replacement  of  the  actual  collection  of  turbules  by  a  collection  of 
simple  turbules  will  not  yield  the  same  result. 


4.  Summary  and  Discussion 


We  describe  the  following  principal  results  in  this  report: 

i;  A  complete  set  of  fluid  equations  involving  the  pressure,  temperature, 
mass  density,  and  flow  velocity  is  presented  for  an  ideal  gas  atmosphere 
for  which  viscosity  and  thermal  conductivity  are  neglected.  These  equa¬ 
tions  include  equations  (8)  or  (9)  as  the  heat  flow  equation. 

ii;  Prom  this  set,  a  complete  set  of  coupled  linear  differential  equations 
is  derived  for  the  acoustic  flow  velocity  and  the  variations  of  acoustic 
pressure,  temperature,  and  mass  density  relative  to  underlying  station¬ 
ary  turbulent  fluid  fields.  Equations  relating  the  turbulent  fields  are 
also  derived. 

iii;  We  show  that  the  coupled  acoustic  differential  equations  yield  Monin's 
wave  equation  for  the  acoustic  pressure  when  only  first-order  terms  in 
the  turbulent  temperature  variation  and  velocity  fields  are  retained. 

iv:  We  show  that  the  length  scale  ag  of  the  scattering  volume  is  just  as 
important  as  the  acoustic  wavelength  A  and  the  length  scales  a  of  the 
turbulence  in  predicting  the  general  behavior  of  the  far-field  acoustic 
scattering  by  turbulence  in  a  first  Born  approximation.  In  particular,  if 
a  <C  Os,  then  the  standard  results  (eq  (48)  and  (49))  for  the  scattering 
amplitudes,  proportional  to  cos^,  result  from  Monin’s  equation  for  all 
a/ A.  If  both  a  >  Us  and  a/ A  >  1,  then  the  scattering  amplitudes  (eq 
(50)  and  (51))  result,  which  are  those  predicted  by  the  short  wavelength 
limit  of  Monin’s  equation,  the  same  as  equations  (48)  and  (49)  but 
without  the  cos  6  factor.  If  A  >  a  >  a^,  other  results  occur. 

v:  For  solenoidal  turbuh  nt  velocity  fields,  we  show  that  both  forward  and 
backward  velocity  scattering  are  essentially  zero  for  all  a/A  if  a  -C  Uj, 
but  nonzero  if  a  >  o^. 

The  first-order  results  of  this  report  for  a  <  have  been  applied  to  scatter¬ 
ing  by  individual  turbules.  In  a  followup  report  [2],  we  will  develop  similar 
results  for  ensembles  of  turbules,  and  a  connection  will  be  made  to  the  Kol¬ 
mogorov  spectrum,  a  Gaussian  spectrum,  and  the  standard  structure  function 
development  for  isotropic  homogeneous  turbulence. 
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In  the  work  reported  here,  we  have  also  found  expressions  for  the  mean 
temperature  and  velocity  differential,  as  well  as  total  scattering  cross  sec¬ 
tions  in  first  Born  approximation  for  a  randomly  oriented  stationary  turbule 
of  a  given  scale  length,  but  arbitrary  morpholog}^  We  have  obtained  and 
compared  expressions  for  these  mean  cross  sections  for  both  Gaussian  and 
exponential  envelope  functions  of  the  turbule  temperature  variation  and  ve¬ 
locity.  We  have  given  a  somewhat  arbitrary  but  general  recipe  for  adjusting 
the  parameters  for  ensembles  of  turbules  with  different  envelope  functions  so 
that  they  have  essentially  the  same  size  and  energy  content.  We  have  shown 
that  the  simplest  possible  model  for  an  individual  turbule  (a  spherically  sym¬ 
metric  structure  spinning  nonuniformly  about  a  fixed  axis)  still  yields  strong 
azimuthal  dependence  in  the  velocity  scattering.  Also,  we  have  shown  that 
we  can  always  find  an  ensemble  of  these  simplest  model  turbules,  with  ran¬ 
dom  orientations  of  the  spin  axes,  that  produces  the  same  ensemble  average 
incoherent  scattering  as  an  ensemble  of  randomly  oriented  turbules  of  ar¬ 
bitrary  morphology.  We  have  also  given  formulas  that  yield  the  equivalent 
simple  turbule  parameters  in  terms  of  any  given  actual  ensemble  averages 
(|r(y)P),  {A,{y)A*{y)). 

In  a  future  report  [2],  we  will  make  use  of  these  results  to  construct  an 
isotropic  homogeneous  ensemble  of  turbules  having  different  scale  lengths 
and  locations,  and  to  derive  the  conditions  for  the  existence  and  the  bounds 
of  the  Kolmogorov  (von  K^man)  spectrum  and  a  Gaussian  spectrum  for 
the  incoherent  scattering  from  such  an  ensemble.  In  later  work,  we  hope  to 
be  able  to  use  the  results  of  this  report  to  treat  the  problem  of  incoherent 
scattering  produced  by  anisotropic  inhomogeneous  ensembles  of  turbules. 
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